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where; 


K12  is  mode  I  fracture  toughness, 


D  is  the  fractal  dimension  . 

Assuming  that  fracture  can  be  modelled  as  a  scaling  fractal,  the 
constant  of  proportionality  of  this  relationship  is  found  to  be  a  product  of 

Young’s  modulus  (E)  and  a  characteristic  length  (ao); 

K|(;=A(D-1)’'2, 

where  A  =  E(ao)^''2_ 

The  constant  A  is  a  family  parameter  which  identifies  a  line  in  the 
toughness-fractal  dimension  plane.  Materials  within  a  given  line  experience  an 
increase  in  tougnners  as  the  fractal  dimension  increases. 

Using  experimentally  determined  values  of  ,  E.  and  D.  a 
cnaracteristic  length,  ag  .  can  be  computea.  This  lengii’,  is  moicative  of  the 

i 

“unit  process"  of  fracture,  and  represents  an  average  step  siae  in  the 
geometry  of  fracture.  For  example,  the  characteristic  length  of  a  polycrystalline 
alumina  (AD90)  is  3  A  ,  suggesting  the  Al-O  bond  rupture  as  characteristic  of 
ihe  fracture.  Alternativety,  a  sine  silicate  glass  ceramic  ( MS^G8«5 )  has  a 
characteristic  length  of  76  A.  This  mdicstos  a  cluster-like  fracture,  suggestive  of 
glass-crystal  mosecular  graypmgs  as  the  unit  process.  Thus,  if  this  mouei  Is 
correct,  v/e  have  a  technique  wnereby  we  car,  relate  the  gccmetiy’  of  the 
fracture  surface  of  brittle  masenais  to  an  atomic  scale  csometry. 
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ABSTRACT 


Fractal  geometry  provides  a  tool  for  the  description  of  irregular  objects. 
While  Euclidean  geometry  allows  for  only  integer  dimensions,  fractal  geometry 
admits  to  the  existence  of  a  dimensional  continuum.  Thus,  geometric  shapes 
can  be  classified  according  to  their  dimensions. 

In  previous  work,  a  relationship  was  shown  to  exist  between  a  materials 
fracture  toughness  and  its  fractal  dimension.  It  was  found  that: 

K|c-(D-1)''2. 

where:  K|c  is  mode  I  fracture  toughness, 

D  is  the  fractal  dimension . 

Assuming  that  fracture  can  be  modelled  as  a  scaling  fractal,  the 
constant  of  proportionality  of  this  relationship  is  found  to  be  a  product  of 

Young’s  modulus  (E)  and  a  cliaracierislic  length  (ao); 

K,c  =  A(D-1)’'2. 
where  A  =  E{3o)'^2, 
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The  constant  A  is  a  family  parameter  which  identifies  a  line  in  the 
toughness-fractal  dimension  plane.  Materials  within  a  given  line  experience  an 
increase  in  toughness  as  the  fractal  dimension  increases. 

Using  experimentally  determined  values  of  K|q  ,  E,  and  D.  a 

characteristic  length,  sq  ,  can  be  computed.  This  length  is  indicative  of  the 

"unit  process"  of  fracture,  and  represents  an  average  step  size  in  the 
geometry  of  fracture.  For  example,  the  characteristic  length  of  a  polycrystalline 
alumina  (AD90)  is  3  A ,  suggesting  the  Al-0  bond  rupture  as  characteristic  of 
the  fracture.  Alternatively,  a  zinc  silicate  glass  ceramic  ( MS498#5 )  has  a 
characteristic  length  of  76  A.  This  indicates  a  cluster-like  fracture,  suggestive  of 
glass-crystal  molecular  groupings  as  the  unit  process.  Thus,  if  this  model  is 
correct,  we  have  a  technique  whereby  we  can  relate  the  geometry  of  the 
fracture  surface  of  brittle  materials  to  an  atomic  scale  geometry. 
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CHAPTER  1 


INTRODUCTION  TO  FRACTALS  AND  FRACTURE 
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Ecaslals 

Fractal  geometry  is  a  non-Euclidean  geometry  that  was  developed, 
popularized  and  applied  by  B,B.  Mandelbrot.^  The  word  fractal  is  derived 
from  the  latin  "fractus"  which  means  fragmented  or  broken.  Mandelbrot  realized 
that  fragmented  geometries  rather  than  Euclidean  geometries  are,  by  far,  the 
most  common  geometries  of  nature.  His  research  has  created  a  framework  for 
the  description  of  systems  that,  up  to  now,  were  too  chaotic  for  geometric 
description. 

In  his  book,  The  Fractal  Geometry  of  Nature.  Mandelbrot  describes  the 
extensive  applications  of  fractal  geometry.  The  concepts  have  been  used  to 

describe  the  geometry  of  clouds,  soot  aggregates,  dielectric  breakdown, 

♦ 

mountainous  terrain,  coastlines,  etc.  This  thesis  uses  the  notions  of  fractal 
geometry  to  describe  the  geometry  of  fracture  surfaces.^ 

Basic  Concents  of  Measure 

To  establish  a  “(eer  for  fractal  geometry  we  can  simply  consider  how  we 
measure  the  length  of  a  line.^  It  is  clear  that  the  measurement  of  a  set  requires 
a  measuring  tool,  or  ruler. 


-SW-V 


Suppose  we  use  a  supply  of  discs  of  a  given  radius,  R,  and  consider  the 
line  of  Figure  1. 


Figui«  1 .  Measuring  the  length  of  a  line  requires  a  ruler 


The  length  of  this  line  is  measured  by  covering  it  completely  with  as  few  discs 
as  possible,  as  in  Figure  2.  The  length,  then,  is  given  by: 


length  ^  (number  of  discs)Xv2R) 
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Now,  suppase  we  have  a  new  supply  of  discs  with  a  different  radius,  say 
(1/2)R.  The  line  length  is  determined  in  the  same  manner,  and  a  new  line 
length  is  computed,  Figure  3,  Note  the  difference  in  the  measured  line  length  in 
Figures  2  and  3.  Although  the  measuring  disc  was  decreased  to  half  of  the 
original  measure,  the  number  required  to  cover  the  line  more  than  doubled. 

The  "tortuous"  nature  of  this  curve  gives  a  length  that  is  scale  dependent.  As 
the  scale  decreases,  the  measured  length  increases. 

The  measured  length  of  a  line  can  be  succinctly  described  by  the 
following  equation,  called  Richardson’s  equation*^ : 

L  =  kE^’D  (1) 

Where:  L  is  the  measured  length, 

E  is  the  measuring  scale, 

D  is  the  dimension  of  the  curve, 
k  is  a  proportionality  constant. 
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If  P  is  equal  to  one.  the  line  is  Euclidean  and  its  measured  length  is  not  a 
function  of  scale.  If  1<D<2,  the  curve  is  said  to  be  fractal  with  a  dimension 
given  by  D.  This  dimension  can  be  thought  of  as  an  indicator  of  the 
"wiggliness"  or  space  filling  nature  of  the  curve. 


If  we  are  given  a  line,  we  can  compute  its  fractal  dimension  from 
equation  (1).  The  line's  length  is  computed  over  a  range  of  scales.  A  log-log 
plot  of  length  vs.  scale  gives  a  straight  line  with  slope  equal  to  1-D,  Figure  4. 


Figure  4.  A  log-log  plot  of  Richardson's  equation 
provides  a  slope  equal  to  i-D. 
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It  should  be  mentioned  that  mathematical  fractals  will  obey  equation  (1) 
for  all  scales.  Thus,  as  the  scale  goes  to  zero,  the  measured  length  becomes 
infinite.  Physical  fractals,  however,  may  have  a  finite  cutoff  at  both  the  large 
and  small  scales.  In  the  physical  world,  the  fluctuation  of  a  line  will  find  its  limit 
in  the  smallest  measureable  featu:-^.  The  sensitivity  of  our  measuring  stick  will 
determine  the  smallest  feature  we  can  measure.  If  fluctuations  occur  below  the 


sensitivity  of  our  instrument,  the  line  will  appear  to  be  Euclidean.  In  this 
instance,  the  upper  bound  of  the  measured  length  of  a  line  is  simply  that  length 
that  corresponds  to  the  finest  possible  scale  of  measurement.  We  could  not 
conclude,  however,  that  the  line  has  exhausted  its  fluctuation,  we  have  simply 
exhausted  our  ability  to  measure  any  fluctuations. 


On  the  contrary,  a  line  may  be  composed  of  indivisible  components. 

After  our  measuring  stick  has  become  smaller  than  the  smallest  of  these 
components,  the  line  would  cease  to  exhibit  scale-dependent  length.  This 
length,  then,  would  represent  a  true  upper  bound.  Such  a  line  could  be  called 
“fractal  over  the  range  it  exhibited  scale-dependent  length.  It  is  not.  however,  a 
fractal  in  the  strict  mathematical  sense. 


Both  of  the  aforementioned  possibilities  would  have  a  Richardson  plot  as 
In  Figure  5.  The  object  would  be  described  as  a  fractal  over  a  bounded  range 
of  scales.  These  cut-off  scales  are  of  special  significance  to  the  particular 
geometry  of  a  specific  physical  phenomenon. 


Figure  5.  Physical  fractals  have  scale  dependence  over 
a  bounded  range  of  scales. 


Fractals  and  Fracture 

Fracture  markings  on  glasses  and  polycryslalline  brittle  materials 
( known  as  mirror,  mist  and  hackle)  are  precursors  to  crack  branching^  and  can 
be  used  to  describe  the  stress  state®  and  characteristics  ot  crack  propagation.^ 
These  markings  have  been  observed  for  over  50  years  and  were  related 


quantitatively  to  the  stress  in  the  1950’s.®  More  recently,  the  repetition  of  these 
features  (multiple  mirrors)  was  observed  and  quantitatively  fGlated  to  stress 
intensity.^  Ravi-Chandar  and  Knauss®  also  recently  noted  that  mist  and  hackle 
are  self-similar;  i.e.,  they  appear  to  be  physically  similar  and  produced  in  the 
same  fashion.  This  observation  follows  previous  detailed  descriptions  of  the 
structure  of  mist  and  hackle.^  ®  These  descriptions  did  not,  hov;ever, 
emphasize  the  self-similar  nature  of  the  features.  These  recent  observations  of 
multiple  mirror,  mist,  hackle,  crack-branching  and  self-similarity  led  to  this 
research. 

A  Basic  Model  of  Fracture 

In  order  to  understand  the  essentials  of  material  fracture,  fracture 
mechanists  employ  a  basic  model  of  fracture,  in  this  model,  atoms  interact  with 
their  nearest  neighbors  via  a  simple  two-body  potential.  The  derivation 
assumes  that  atoms  lie  in  a  plane,  that  they  interact  only  with  their  nearest 
neighbor  and  that  bonds  break  sequentially.  That  is.  fracture  is  localized  to 
only  two  atoms  at  a  time. 


The  potential  energy  between  atoms  in  a  crystal  structure  can  be 
approximated  by  a  function  of  the  form"*  ; 

U  = -Qi/ r^  +  Qg/ rV,  (2) 

where  -Ql /  r^  is  an  attractive  term, 

Q2/  rV  is  a  repulsive  term. 

For  ionic  crystals,  the  exponent  in  the  attractive  term  (x)  is  close  to  unity  while 
the  repulsive  exponent  (y)  is  between  6  and  12.  Such  a  function  would  have 

the  form  shown  in  Figure  6,  where  dg  is  the  equilibrium  spacing  between 

neighboring  atoms.  The  depth  of  the  potential  minimum  is  related  to  fracture 
energy,  while  the  curvature  of  the  potential  minimum  is  an  indication  of  the 
modulus.  The  derivative  of  this  potential  function  gives  the  force  function 
between  neighboring  atoms,  as  shown  in  Figure  7. 


Figure  6.  Potential  energy  between  neighboring  atoms 


stress 

Theoretical  cohesive  strength 

. . 

Figure  7.  Force  function  between  neighboring  atoms. 


When  a  material  is  stressed  beyotiid  a  certain  limit,  irreversible  changes 
wit’  occur.  If  loading  is  continued,  the  material  can  be  fractured.  As  a  crack 
propagates  through  a  specimen,  the  potential  energy  of  the  specimen  is 
reduced. 
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Energy  is  released  via  the  creation  of  new  surface,  the  emission  of  light, 
acoustic  energy,  etc.  The  energy  required  to  propagate  a  crack  is  called  the 

fracture  energy,  y^. 

Crack  propagation  through  a  microstructure  is  a  geometrically 
complicated  statistical  process  that  contains  atomic  scale  as  well  as 
macroscopic  scale  geometries. ^  ^  Interaction  of  the  crack  front  with  the 
microslructure  will  affect  the  geometry  and.  hence,  the  surface  area  generated. 

A  material's  resistance  to  crack  propagation  is  called  fracture  toughness. 
Fracture  toughness  may  be  characterized  by  the  fracture  energy,  ,  the  strain 

energy  release  rate,  »  2y|).  or  the  critical  stress  intensity  factor 

(  Kq  »  v'  2Ey|  ).  A  toughening  mechanism,  then,  is  a  microstructural  property 

that  affects  a  material's  fracture  toughness.  One  particular  class  of  tougirening 
mechanisms  relies  on  influencing  the  direction  of  crack  propagation 
(e.g,,  crack  deflection  An  increase  in  toughness  can  be  achieved  by 
making  the  crack  path  irreguiar.^^  Thus,  the  geometry  of  the  Iraciure  surface 
plays  a  role  in  the  material's  toughness. 


The  importance  of  fracture  energy  is  made  obvious  in  the  following 
simplified  derivation  assuming  plane  fracture  in  a  perfect  crystal.^  ^  .15,16 
A  crack  is  assumed  to  reside  between  two  planes  of  atoms  with  interplanar 

spacing dg,  Figures. 


Raihor  than  using  the  derivative  of  eQuaUon  (2).  the  strep's  resulting  from  the 
strain  of  the  inlerpianar  bonding  is  approximated  by  a  sinusoidal  function. 
Figure  9. 
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Figure  9.  Sinusoidal  approximation  of  the  force  between 
neighboring  atoms. 


At  the  equilibrium  spacing,  dQ,  the  stress  is  zero,  rising  to  a  maximum 
and  dropping  to  zero  at  a  distance  Jq+h,  the  distance  of  bond  rupture. 

Therefore,  o=  Oq  sin  {7t(x*do)/a),  (3) 

where;  x  is  the  separation  beyond  equilibrium, 

a  is  the  failure  separation. 

For  small  displacements,  x,  the  sin  is  equal  to  its  argument. 


Therefore,  o  =  Oq  it  (x-do)/a  , 


(4) 
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By  definition,  the  strain  is; 


e  =  x/do , 


and,  for  small  displacements, 


a  =  Ee  =  Ex/dg , 


so  that 


da/dx  =  E/dg . 


Comparing  equations  (5)  and  (8)  provides; 


Og  =  Ea/jidg . 


'5 


*  ^.>^1 

* 

Sr 


Thus,  an  estimate  of  the  theoretical  cohesive  strength  is  obtained  from  known 
values  of  E,  a .  and  dg.  Now,  in  fracture,  the  work  of  separation  must  at  least 

be  equal  to  the  energy  required  to  create  two  new  surfaces,  ^  .15.16 


The  work  per  unit  of  area  generated  by  fracture  is; 

work/unit  area  =  J  OQSinjt(x-do)/a  dx  ==  ,  (10) 

so  that,  ZcQdJn  =  2y^,  (11) 

but,  from  equation  (7); 

2Ea2/]i2dQ  =  273  ,  (12) 

SO  YgsEa2/it2do.  (13) 

Equation  (13)  shows  that  the  surface  energy  tor  plane  fracture,  Yg ,  is  a 
function  Of  modulus.  E.  interplanar  spacing,  dg.  and  the  displacement  to  bond 

rupture,  a.  For  single  crystals,  the  surface  energy,  ,  is  assumed  to  be  nearly 

equal  to  the  fracture  enrgy,  Yj .  Gilman  1®  assumed  that  the  constant,  a,  was 
equal  to  the  radius  of  atoms  within  a  given  plane. 


For  single  crystals,  his  results  were  uniformly  high,  yet  of  the  same  order  of 
magnitude  as  experimentally  obtained  values.  Polycrystalline  fracture 
energies,  however,  are  approximately  an  order  of  magnitude  higher  than  their 
single  crystal  counterparts.  ^  ^  Part  of  this  difference  is  due  to  an  increase  in 
the  roughness  of  the  fracture  surface. 

The  model  of  fracture  presented  here  takes  no  account  of  fracture 
surface  roughness.  It  is  assumed  that  fracture  creates  two  perfectly  plane 
fracture  surfaces.  Fractal  geometry  can  accommodate  surface  roughness  by 
allowing  for  irregular  paths.  These  paths  could  be  constructed  from  straight  line 
segments  imbedded  in  an  orderly  geometric  matrix. 

Imagine  a  two-dimensional  array  of  atoms  as  shown  in  Figure  10. 
Fracture  could  progress  from  region  A  to  region  B  through  a  number  of  paths. 
Some  of  these  paths  are  shown  in  the  Figure  ( a-a,  b-b,  and  c-c).  The  basic 
model  would  assume  only  a  straight  line,  whereas  a  fractal  model  could 
describe  Other  more  complicated  paths.  Thus,  surface  rougtmess  can  be 
accommodated  without  a  radical  departure  from  the  basic  model. 


•  « *  o  •fc* 

4.  *» « *  r 


•A  •a  V-  a  /  -  f  J* 


HE 

H 

S 

jip 

3 

H 

1 

jS 

3 

mm 

1 

H 

gs 

1 

m 

Figure  1 0.  Schematic  of  possible  crack  paths  in  an  idealized  lattice. 
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Fr,aclure_Enerqy  and  Fracture  Toughness 


Equation  (9)  provides  an  estimate  of  the  theoretical  cohesive  strength  of 
a  material.  Brittle  ceramics  exhibit  fracture  strengths  that  are  generally  orders 
of  magnitude  lower  than  this  estimate.  This  difference  results  from  the  stress 
concentrating  effects  of  flaws. 


Griffith^  2  able  to  calculate  the  stress  at  which  a  crack  would 


propagate.  He  postulated  that  crack  extension  would  occur  only  if  the  creation 
of  new  surface  acts  to  reduce  the  potential  energy  of  the  system.  He  found  that, 
under  conditions  of  plane  stress; 


(14) 


Where:  Oj  is  the  fracture  strength, 


E  is  Young’s  modulus, 


is  the  fracture  energy, 


C  is  the  flaw  size. 
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Inglis  made  computations  of  the  stress  concentration  created  by  flaws  of 
various  geometries.  It  was  found  that  cracks  in  a  material  could  be 
conveniently  modelled  as  elliptical,  and  their  stress  concentration  computed. 
The  stresses  near  the  tip  of  such  elliptical  cracks  have  a  precise  mathematical 

description  that  depends  upon  a  parameter  called  stress  intensity,  K| 

( The  I  denotes  tensile  loading).  For  a  material  to  fail,  the  far  field  stress  must 
be  concentrated  to  the  level  of  theoretical  cohesive  stress.  Thus,  stress 
intensity  reaches  a  critical  value  called  the  fracture  toughness,  or  critical  stress 

intensity  factor,  K|q  .  This  critical  value  is  found  to  depend  on  the  flaw  size  and 
the  far  field  stress: 

K|c=Y(j^a/C.  (15) 

Where:  oj  is  the  far  field  stress, 

C  is  the  Haw  size. 

Y  is  a  factor  that  depends  on  the  geometry  of 


loading  and  crack  configuration. 
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Comparing  equations  (14)  and  (15)  provides  a  relationship  between  K|q  and  y^. 
Under  conditions  of  plane  stress; 

K,o  =  (2E7,)1/2  (16) 


Fracture  Profiles  as  Fractal  Curves 

In  many  cases,  a  fracture  surface  presents  a  rather  complicated 
geometric  structure.  This  complex  geometric  shape  is  not  amenable  to 
description  by  Euclidean  geometry,  thus,  fracture  surface  energies  have  not 
been  accurately  modeled  as  functions  of  this  geometry.  Fractals,  however, 
provide  a  tool  for  the  description  of  such  surfaces.  At  the  very  least,  they  can  be 
used  to  quantify  the  surface  roughness.  A  fracture  profile  presents  a  curve 
similar  to  that  of  Figure  1 1 .  The  "complexity"  of  this  profile  can  be  categorized 
by  its  fractal  dimension.  As  discussed  earlier.  Richardson's  equation,  where 
the  profile  length  is  measured  for  various  scales,  provides  a  method  lor 
determining  the  fractal  dimension.  Of  course,  a  fracture  surface  is,  at  the 
very  least,  a  two-dimensional  object,  but  as  a  first  step  in  categorization,  it  is 
worthwhile  to  compute  the  dimension  of  the  profile,  Figure  1 1 . 


Fractured  sample 


Profile  view 


I _ I 

Figure  1 1 .  Fracture  profiles  can  be  measured  according  to 
Richardson. 
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It  is  quite  possible  that  there  will  be  a  physical  limit  to  the  scale 
dependence  of  length.  Features  on  the  fracture  surface  may  vary  only  down  to 
a  finite  cut-off  size.  This  may  be  on  the  grain  size,  molecuiar,  or  atomic  scale. 
But,  this,  in  itself,  is  useful  information. 

Scaling  Fractals 

There  is  a  particular  subclass  of  fractal  geometry  which  is  comprised  of 
"scaling  fractals".^  This  is  what  is  generally  referred  to  as  a  fractal. 
Regardless  of  the  nature  of  a  curve.  Richardson’s  equation  can  be  used  to 
compute  a  fractal  dimension.  Scaling  fractals,  however,  have  two  very 
important  properties:  scale  invariance,  and  seif-similarity.  Scale  invariance 
refers  to  features,  geometrically  identical,  appearing  on  all  scales  of 
observation.  Self-similarity  means  that  a  small  feature  of  the  object  can  be 
scaled  to  precisely  match  a  large  feature.  These  concepts  are  most  easily 
understood  by  obsen/ing  the  construction  of  a  scaling  fractal. 

Consider  a  box.  as  shown  in  Figure  12a^'^'^^,  whose  sides  are  of  unit 
length.  Next  to  this  box  is  a  shape  called  a  generator,  composed  of  tine 
segments  scaled  to  1/4  of  the  length  of  a  side  of  the  box.  Now.  replace  each 
side  of  the  square  with  this  “scaled  ”  generator. 
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The  result  is  shown  in  Figure  12b.  Again,  take  the  generator  and  scale  it  down 
to  the  length  of  each  straight  line  segment  of  this  new  object.  The  generator  is 
now  composed  of  line  segments  of  length  (1/4)2,  Replace  each  of  the  straight 
line  segments  of  Figure  1 2b  with  the  scaled  generator.  A  portion  of  the 
resulting  shape  is  shown  in  Figure  12c.  This  process  is  continued,  ad  infinitum, 
generating  an  object  which  is  said  to  be  scale  invariant  and  self-similar. 


Each  side  ot  the  square 
is  replaced  by  the  scaled 
generator. 


Figure  12b 


Figure  12c 


Figure  12.  The  stepwise  construction  ot  a  scaling  fractal. 
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A  scaling  fractal  looks  geometrically  the  same  everywhere  and  on  all 
scales.  Clearly,  we  have  generated  an  object  whose  perimeter  will  be  a 
function  of  scale  (as  in  equation  1).  The  fractal  dimension  of  this  object, 
however,  can  be  computed  in  another,  quite  different  fashion.”*  Scaling 
fractals  obey  the  following; 


NrD  =  i. 


(17) 


Where:  N  is  the  number  of  elements  in  the  generator, 

r  is  the  scale  factor  of  an  element. 

D  is  the  similarity  or  fractal  dimension. 

For  the  curve  of  Figures  12,  called  the  quadratic  Koch  curve,  N  =*  8,  r  a  i/4. 
Thus. 

D  logs /log4  =  3/2.  (18) 


Figure  13  contains  three  photographs  of  a  glass  fracture  surface.  Such 
fractures  contain  readiSy  identifiable  regions  called  mirror,  mist,  and  hackle.  A 
glance  at  the  photographs  will  demonstrate  the  similarity  of  the  mist  and  hackle 


regions. 


Figure  13b  Is  a  photo  of  the  mist-hackle  boundary,  while  13a  is  a  view 
of  the  mirror-mist  boundary.  Figure  13c  is  a  magnified  view  of  a  region  in  13a. 
The  identification  of  these  boundaries  is  clearly  affected  by  magnification. 

Given  such  scale  dependence,  it  seems  possible  that  a  geometry  that 
describes  one  region  may  simply  be  a  scaled  version  of  a  geometry  that 
describes  another  region.  This  is  precisely  the  notion  of  scaling  fractals  and 
fractal  geometry. 

It  is  certainly  true  that  fracture  profiles  can  be  measured  according  to 
Richardson's  equation  and.  if  such  is  the  case,  assigned  fractal  dimensions. 

The  proposal  of  this  thesis  is  that  fracture  surfaces  can  be  accurately  modelled 
as  scaling  fractals.  Thus,  we  can  imagine  a  scheme  of  generation  and  a 
generator  that  is  repeated  and  scaled  as  fracture  progresses.  This  generator 
may  be  an  atomic  scale  process  that  cascades,  through  the  scale  invariance  of 
fractal  geometry,  to  the  macroscopic  features  of  a  fracture  surface.  Given  D.  we 
can  use  equation  14  together  with  knowledge  of  the  crystal  structure  of  a 
material  to  determine  possible  values  for  r  and  N.  Modelling  fracture  in  this 
fashion  allows  for  geometric  interpretation  of  measurements  on  the 
macroscopic  scale,  to  a  description  of  geometry  on  the  atomic  scale. 
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Research  has  shown  a  relationship  between  fracture  toughness  (i.e., 
resistance  to  crack  growth)  and  fractal  dimension.'*  >24  At  first 

glance,  this  may  not  seem  surprising.  After  all,  fracture  toughness  is,  in  some 
ways,  related  to  surface  roughness,  and  fractal  dimension  provides  a  measure 
of  this  roughness.  The  relationship  is  not  simple,  however,  and  contains  some 
interesting  results. 


CHAPTER  2 

EXPERIMENTAL  TECHNIQUE 


In  this  study,  two  independent  techniques  were  used  to  determine  the 
fractal  dimension  of  fractured  samples:  slit-island  analysis^  and 
Richardson  plots.  Passoja^^  used  a  third  technique,  Fourier  transform 
analysis,  on  a  zinc  silicate  glass  ceramic  and  obtained  values  within  2%  of 
those  obtained  by  slit-island  analysis  and  Richardson  plots. 


Slit-Island  Analysis 

Fractured  samples,  obtained  from  previous  studies, were  carefully 
cleaned  and  coated  with  nickel.  The  nickel  coating  performs  two  functions:  It 
provides  good  contrast  during  polishing,  and  it  helps  to  hold  the  fracture 
surface  together.  The  samples  are  then  potted  in  epoxy  and  polished  parallel 
to  the  fracture  surface,^  ^<22  Figure  1 4.  As  the  fracture  surface  is 
encountered,  a  section  of  the  fracture  surface  appears  in  the  polishing  plane. 
These  sections  appear  as  islands  in  the  polishing  plane.  As  polishing 
proceeds  these  islands  begin  to  grow.  The  perimeter  of  the  islands  presents 
line  that  can  be  measured  according  to  Richardson’s  equation. 


Figure  14.  Fractured  samples  are  encapsulated  in  epoxy  and  polished  parallel 
to  the  fracture  plane.  Islands  emerge  in  the  polishing  plane. 

An  adaptation  of  Richardson’s  equation  reveals  a  relationship  between 
the  area  and  perimeter  of  an  island: 

A~p2/D  (19) 

where;  A  is  the  island  area, 

P  is  the  island  peiimeter, 

D  is  the  dimension  of  the  perimeter. 


At  the  first  emergence  of  an  island,  a  Polaroid  photograph  is  taken  at 
magnifications  ranging  from  10X  to  400X.  Polishing  proceeds  in  stages,  each 
stage  foiiowed  by  a  photograph.  Once  a  magnification  is  chosen,  that  same 
magnification  is  used  for  all  subsequent  photos.  As  polishing  progresses,  the 
isiands  grow  and  sometimes  merge.  A  representative  sequence  of  photos, 
taken  at  10X  magnification  is  shown  in  Figure  15.  These  photos  are  of  a  zinc 
silicate  glass  ceramic  (MS-508). 

The  "roughness"  of  the  perimeter  is  an  indication  of  the  roughness  of 
the  f  racture  surface.  A  sequence  of  approximately  20  photographs  is  taken  to 
document  island  growth.  These  provide  data  fci  a  log-log  plot  of  area  vs. 
perimeter  (equation  1 9),  from  which  a  slope  of  2/D  is  obtained.  A 
representative  graph  is  shown  in  Figure  1 6.  Data  and  graphs  for  the  materials 
tested  are  contained  in  following  sections. 


Photo  #4 


Photo  #7 


Photo  #8 


Figure  15.  Sample  photographs  from  zinc  silicate  MS508 
slit-island  analysis. 


Figure  16.  Fractal  dimension  can  be  obtained  from 
area-perimeter  growth  of  slit-islands. 


Occasionally,  an  island  will  encounter  the  machined  edge  of  a  sample, 
(e.g..  Figure  1 5).  The  island  will  have  a  perimeter  that  is  composed  of  two 
sections:  A  line  representative  of  the  fracture  surface,  and  a  line  representative 
of  the  machined  edge.  The  resulting  fractal  dimension  would  be  a  weighted 
average  of  these  two  portions.  Yet,  our  interest  is  only  in  that  portion  that 
represents  the  fracture  surface.  Such  edge  encounters  are  simply  folded  out  of 
existence.  Figure  1 7  explains  graphically  how  the  island  can  be  folded  across 
the  edge  to  create  an  island  with  twice  the  original  area  and  twice  the  "fracture 
perimeter."  Thus,  we  are  left  with  an  island  that  is  unbiased  by  the  dimension 
of  the  sample  edge. 


Island  areas  and  perimeters  are  computed  from  the  polaroid 
photographs  using  a  Lemont  Scientific  image  analysis  system."*  ^  Different 
grey  levels  of  the  photographs  are  assigned  different  colors.  The  system 
software  is  then  able  to  compute  area  and  perimeter  of  the  different  colored 
regions. 

In  brief,  the  system  places  a  grid  with  one  micron  spacing  atop  the 
selected  image.  Points  of  intersection  of  the  grid  and  island  perimeter  are  then 
assigned  Cartesian  coordinates.  Straight  lines  are  drawn  between  these 
points  and  simple  trigonometry  allows  for  the  computation  of  the  length  of  these 
lines.  The  island  perimeter  is  taken  to  be  the  sum  of  these  straight  line 
segments.  Island  area  is  equal  to  the  area  of  all  totally  enclosed  rectangles 
plus  the  triangular  regions  along  the  edges.  Figure  1 8  demonstrates  this 
procedure. 


1  micron 


t 


Cartesian  coordinates  are  assigned  to  the 
intersection  points  of  the  grid  and  island 


Straight  lines  connect  the  intersection  points. 
The  area  and  perimeter  are  computed  from  this 
approximation. 


Figure  18.  A  grid  is  placed  atop  the  island  of  interest.  Intersection  points  of 
the  grid  and  perimeter  are  used  to  compute  area  and  perimeter. 
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The  area  and  perimeter  information  is  used  to  plot  a  graph  for  each 
island  that  appears  on  the  fracture  surface.  Log-log  plots,  as  in  Figure  16,  are 
used  to  extract  the  fractal  dimension. 

As  a  check  on  the  enliro  procedure,  an  object  of  known  dimension  was 
analyzed  with  the  slit-island  technique.  The  chosen  object  was  an  alumina 
sphere.  Any  cutting  plane  through  a  sphere  reveals  a  cl.xle,  the  dimension  of 
which  is  equal  to  one.  The  sphere  was  nickel  coated,  encapsulated  in  epoxy 
and  polished.  Sections  through  the  sphere  were  photographed  at  1 0X,  Figure 
19,  and  image  analyzed.  The  results  are  shown  in  Figure  20.  The  computed 
dimension  was  equal  to  .98,  a  difference  of  only  2  percent  from  the  expected 
value  of  a  perfect  sphere,  0*1 .0. 
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Photo  #1 


Photo  »^10 


Figure  19*  S'  ..pie  photographs  from  slit-island  analysis 
of  an  alumina  sphere. 


y  =  '  3.0657  +  2.0338X  R  =  0.99 


togP 


This  graph  shows  the  experimental  data  obtained 
through  slit-island  analysis  ot  a  Euclidean  obiect. 
Such  an  object  should  have  a  dimension  equal  to 
one.  This  graph  results  in  a  dimension  oi  .98 


Figure  20,  Data  from  the  Euclidean  test  specimen. 
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Richardson  Plot 

Richardson  plots,  as  mentioned  in  the  introduction,  detail  the  change  in 
measured  length  of  a  line  as  a  'unction  of  scale.  Therefore,  construction  of  a 
Richardson  plot  for  a  fracture  surface  requires  a  line  that  is  representative  of  the 
fracture  surface  and  a  range  of  scales  for  measuring  that  line.  As  in  slit-island 
analysis,  fractured  samples  are  coated  with  nickel  and  encapsulated  in  epoxy. 
The  sample  is  then  polished  either  parallel  or  perpendicular  to  the  fracture 
plane.  If  polished  parallel,  we  obtain  an  island  whose  perimeter  is 
representative  of  the  fracture,  if  polished  perpendicular,  we  obtain  a  fracture 
profile.  In  either  case  a  photo  montage  is  constructed  from  polaroid 
photographs  at  400X  magnification.  This  montage  is  then  measured  with 
dividers  set  to  2,  1/2,  1/4,  and  1/8  mm  openings,  as  depicted  in  Figure  21. 
The  length  of  the  profile  (or  perimeter)  is  tabulated  for  each  divider  setting.  In 
this  fashion,  profile  length  is  computed  as  a  function  of  scale.  A  log-log  plot  of 
length  vs.  scale  (Figure  4)  gives  a  straight  line  with  slope  equal  to  1-D. 


Figure  21 .  A  montage  of  the  profile  is  measured  with  dividers 


The  scale  (or  divider  opening)  is  a  measure  of  discernibility.  The 
smallest  observable  feature  is  that  of  the  scale.  As  the  scale  becomes  finer, 
we  observe  greater  ar\d  greater  detail.  If  the  scale  is  larger  than  the  largest 
features,  the  observation  is  insensitive  to  those  features.  Thus,  a  curve  will 
begin  to  look  fractal  only  after  the  scale  becomes  smaller  than  such  features. 
Figure  22  shows  a  sample  montage  for  a  section  of  island  perimeter  of  an 
Ocala  chert  ( chert  U5). 


CHAPTER  3 


RESULTS  AND  DISCUSSION  • 


Mecholsky,  Feinberg  and  Passoja^^  determined  a  relationship 
between  tracture  toughness  (Kj^)  and  fractal  dimension  (D).  They  found  that; 

K|o-(D-1)'/2  (20) 

Thus,  as  fractal  dimension  increases,  fracture  toughness  increases. 

(The  term  (D-1 )  will,  henceforth,  be  referred  to  as  D* .  D* ,  then,  is  the 
fractional  part  of  the  fractal  dimension).  They  obtained  the  fractal  dimension  of 
a  number  of  zinc  silicate  glass  ceramics  and  aluminas.  Table  I,  Figure  23. 

The  data  conformed  to  a  straight  line  in  the  log-log  plane  with  a  correlation  of 
.92. 

A  number  of  new  materials  were  chosen  for  testing.  A  complete 
summai7  of  these  materials,  the  technique  used  for  the  determination  of  their 
fractal  dimension,  and  the  resulting  fractal  dimension  is  contained  in  Table  II. 

A  complete  listing  of  the  individual  data  tables  and  graphs  follows  this  section. 

The  Ocala  cherts,  commonly  referred  to  as  flint,  are  similar  in  structure 
to  glass  ceramics  (crysiallites  imbedded  in  an  amorphous  matrix).  These 
cherts  represent  a  series  of  heat  treatments  that  alter  the  fracture  toughness 
and  fracture  surface  tci.'ography. 


U5  chert  was  untreated  and  exhibits  the  highest  toughness.  Numbers  three, 
four  and  F5  chert  were  heat  treated  at  3000,  4000,  and  SOOC^S  respectively. 


It  is  common  to  assume  that  fracture  in  single  crystals  will  occur  along  a 
single  plane.  In  general,  loading  conditions  and  inherent  flaws  will  alter  the 
planar  propagation  of  a  crack  through  such  crystals.  It  is  not  certain,  then,  that 
the  dimension  of  the  fracture  surfaces  of  single  crystals  will  always  be  equal  to 
one  (plane  fracture).  Single  crystals  of  Calcium  Fluoride  and  Spinel  were 


chosen  to  demonstrate  this  possibility.  Large  grain  (-SOOiim)  CdTe  and  ZnSe 
were  dimensioned  as  a  first  step  in  determining  grain  size  effects  on  the  fractal 
dimension.  Pyroceram  9606,  a  magnesium-aluminosilicate  material,  and 
IVISSOS  Zinc  Silicate  are  glass  ceramics. 
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TABLE  I 

Data  from  Mecholsky,  Feinberg  and  Passoja'*^'^^ 


Material 

Era.c.tur.e.T.QUfl!ing.s.s. 

K|c  { MPa 

Fractal  Dimension 

MS  498  6b  #4 

1.6 

1.05 

MS  50012b  #5 

1.8 

•  1.09 

MS  500  #5 

2.2 

1.11 

MS  498  #5 

2.0 

1.16 

LAS  glass  ceramic 

2.7 

1.18 

Monsanto  AI2O3 

3.5 

1.21 

GE  AI2O3 

3.9 

1.23 

Lucalox 

4.0 

1.31 

»  M  •  «  •  >  • 
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D‘ 


Figure  23.  A  graph  of  the  data  of  Table  I  shows;  K|q~D*''^2 


TABLE  II 


Data  Summary 


Mataiiial 


CdTe 

.7 

1.2 

S.l. 

ZnSe 

.9 

1.3 

S.l. 

PZT 

1.5 

1.23 

R 

MS  498 

1.6 

1.05 

S.l. 

MS  500#4 

1.8 

1.09 

S.l. 

MS  500#5 

2.2 

1.11 

S.l. 

MS  498#5 

2.0 

1.16 

S.l. 

Alumina  #3 

3.5 

1.21 

S.l. 

Alumina#5  GE 

3.9 

1.23 

S.l. 

Alumina#3  GE 

3.9 

1.06 

S.l. 

Pyroceram 

2.4 

1.17 

s.l. 

Spinel(s.c.) 

1.2 

1.09 

s.l. 

CaF2  (s.c.) 

.3 

1.07 

s.l. 

Poly  Spinel 

2.1 

1.13 

s.l. 

U5  chert 

1.55 

1.32 

R.S. 

#3  chert 

1.46 

1.26 

R 

#4  chert 

1.25 

1.24 

R 

F5  chert 

1.05 

1.17 

R.S. 

*S.I.  refers  to  the  slit-island  measurement  technique,  R  relers  to  Richardson's  technique 
s.c.  Is  a  single  crystal  sample 


Analysis  of  these  new  materials,  at  first  glance,  seemed  to  indicate  no 
definite  relationship  between  fractal  dimension  and  fracture  toughness. 

Indeed,  as  shown  in  Figure  24,  the  data  is  thoroughly  scattered  in  the  KjQ-D* 

plane.  If,  however,  we  assume  that  the  relationship  of  Equation  20  is  correct, 
we  can  order  the  data. 

PROPOSITION:  There  exists  a  family  of  lines  in  the  log-log  K|q-D*  plane  all  of 
slope  1/2. 

A  graphical  representation  of  this  proposition  is  shown  in  Figure  25. 
Though  we  have  no  mathematical  proof  of  this  proposition,  we  may  be  able  to 
demonstrate  "reasonableness." 

The  existence  of  these  curves  implies  a  relationship  between  the  points 
within  any  given  line.  Such  points  constitute  a  family  and  will  be  characterized 
by  a  family  parameter.  Any  one  of  these  lines  will  have  the  following  functional 
form; 

K|c  =  AD*‘'/2  (21) 


where; 


A  is  a  family  parameter. 
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The  family  parameter,  A,  identifies  a  line  in  the  Kjq-D*  plane.  Within  a  family, 

an  increase  in  fractal  dimension  corresponds  to  an  increase  in  fracture 
toughness. 

A  vertical  line  in  the  K|q-D*  plane  will  intersect  a  number  of  family  lines. 

Hence,  materials  which  exhibit  the  same  fractal  dimension  do  not  necessarily 
have  the  same  toughness.  Similarly,,  materials  of  equal  toughness  do  not 
necessarily  have  the  same  fractal  dimension.  Thus,  the  functional  relationship 
of  Equation  (21 )  is  useful  only  in  the  comparison  of  materials  within  the  same 
family. 


Figure  24.  The  data  are  apparently  scattered  in  the  K|q*D*  plane 


55 


The  Toughness-Fractal  Dimension  Relationship 

As  noted  previously,  toughness  is  proportional  to  fractal  dimension; 


Klc-D' 


1/2 


(22) 


The  parameter  which  will  identify  a  family  in  this  space  is  the  constant  of 
proportionality,  A; 


K|c  =  AD*^/2  (23) 

Since  0*  is  dimensionless,  a  dimensional  analysis  of  this  equation  requires 
that  A  have  the  dimensions  of  toughness.  Suppose  that  this  constant  can  be 
further  reduced  to  one  involving  Young's  modulus,  an  indicator  of  bond 
strength,  and  some  characteristic  length  parameter.  The  rationale  for  this 
supposition  stems  from  the  special  importance  attributed  to  scaling  fractals. 
Scaling  iraclals  are  constructed  from  a  unit  process;  i.e..  a  generator-shape 
which  is  scaled  according  to  a  scheme  for  generation.  Thus,  the  geomehY  of 
fracture,  if  it  can  be  accurately  n^odeied  by  scaling  fractals,  will  have  such  a  unit 


uM 


m 


ik' 


process. 
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We  propose  that  the  fracture  topography,  in  its  atomic  form,  is  recorded 
on  the  fracture  surface  through  the  scaiing  cascade  of  fractal  geometry.  There 
wili,  then,  be  an  atomic  scale  generator  with  a  characteristic  length.  Suppose 
tiie  constant  A  is  a  product  of  this  characteristic  length  and  Young's  modulus, 

sothat  A  =  E(aQ)'*'^2. 

Then; 

K|o=  E  (aoD-)1«,  (24) 

“here:  l5  (1,0  fracture  toughness, 

E  is  Young's  modulus, 

D*  is  the  fractal  dimension, 

ag  is  the  characteristic  length. 

The  plane  fracture  model  of  Chapter  1  demonstrated  the  importance  of 

thermodynamic  surface  energy,  ,  in  fracture.  This  derivation  assumed  ideal 

plane  fracture.  The  real  world,  however,  is  far  from  ideal.  Other  processes 
will  absorb  energy,  the  most  obvious  of  which  is  the  non^planar  propagation  of 

a  crack.  A  new  quantity,  called  the  fracture  energy,  y^ ,  is  used  as  an  effective 
fracture  energy. 
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This  term  includes  processes  that  consume  energy  by  the  formation  of  surface, 
zones  of  plastic  deformation,  acoustic  and  photo  emission,  etc. 

Under  conditions  of  plare  stress,  fracture  toughness  is  related  to  fracture 
energy; 


K|e  =  (2EY,)1'2, 


so  that  K|(;2  =  2EV(. 


Vyhere;  Yf  'S  the  fracture  energy. 


Now,  from  equation  23: 


(25) 


(26) 


K,c2  =  ESagD-. 


(27) 


k  *  4i**'***«^«'*«.*k'’  . 

*4  *  ^ 
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Inserting  equation  27  uuo  equation  26  gives: 

(Y{/E)  =  1/2aoD*.  (28) 

Equation  (28)  can  now  be  used  to  compute  the  value  of  bq.  Using 

experimentally  obtained  values  of  yj  and  E,  and  the  experimentally  determined 

value  of  D,  a  plot  of  (Y|/E)  vs  (agD)  is  constructed.  Figure  26.  This  plot  is’ 
derived  to  have  a  slope  of  1/2.  The  materials  that  were  scattered  throughout 
the  K|q-D*  plane  of  Figure  25  are  now  situated  along  a  single  line. 

The  macroscopic  fracture  energy,  yj,  is  normalized  to  Young’s  modulus. 

E.  Young's  modulus  can  be  conceptualized  as  a  measure  of  atomic  bond 
strength.  It  represents  a  weighted  average  of  all  the  possible  bond  strengths  of 

the  material.  Thus,  as  Yf/E  increases,  the  fracture  energy  is  increasing  with 

respect  to  E.  Along  the  horizontal  axis,  ag  represents  an  atomic  scale  length 
while  D*  represents  a  measure  of  surface  roughness.  Moving  outward  along 
this  axis  demonstrates  another  interplay  between  microscopic  (ag)  and 


macroscopic  (O’)  properties. 


More  importantly,  the  vertical  axis  represents  energy,  the  horizontal  axis 
represents  geometry.  This  implies  an  intimate  connection  between  energy  and 
geometry.  .Equation  (28)  provides  a  clear  statement  of  this  connection.  UvSing 
fractal  geometry,  we  may  be  better  able  to  compute  the  actual  fracture  energies 
from  fractal  theoretic  models. 

Table  III  documents  the  calculated  values  for  the  characteristic  length, 
ag,  for  a  number  of  materials.  Notice  that  the  single  crystals  and  large  grain 

potycrystals  (ZnSe,  CdTe)  have  characteristic  lengths  of  a  few  angstroms. 
These  values  are  of  the  order  of  lattice  parameters  and  suggest  that  the  unit 
process  of  fracture  for  these  materials  is  on  the  order  of  atomic  bond  breaking. 
The  zinc  silicate  glass  ceramics  show  much  larger  characteristic  lengths.  This 
suggests  that  fracture  of  these  materials  is  a  cluster*like  or  molecular  process. 
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TABLE  III 


Summary  of  tho  Calculated  Values  of  Characteristic  Length  gq 


EXQPa\  .  J3.  SqJM 


Alumina 

AD90 

11 

390 

1.21 

3 

AD9S9 

19 

406 

1.31 

3 

Lucalox 

26 

305 

1.31 

5 

Zinc  Silicates 

MS  498#5 

27 

90 

1.07 

76 

MS  500#4 

27 

89 

1.12 

53 

Pyroceram 

25 

120 

1.17 

20 

CdTe 

1; 

40 

1.20 

3 

ZnSe 

5.5 

69 

1.30 

4 

CaF2{s.c.)* 

0.5 

114 

1.07 

1 

Spinel{s.c.) 

3 

240 

1.09 

3 

s.c.  refers  to  single  crystal 


The  determination  of  this  characteristic  length,  sq,  is  of  startiing 

significance.  This  length  came  about  from  a  number  of  important  assumptions. 
First,  it  is  assumed  that  fracture  is  a  fractal  process.  Further,  it  is  assumed  that 
this  fractal  process  is  scale  Invariant  and  self-similar.  From  this  assumption, 
we  propose  that  fracture  on  the  macroscopic  scale  is  composed  of  the  scaied 
insertion  of  a  unit  process  on  the  atomic  scale.  In  other  words,  there  is  some 
small  scale  event  with  a  definite  geometric  shape.  This  shape  is  repeated 
everywhere  along  the  fracture  and  reveals  its  signature  from  the  atomic  to  the 
macroscopic  scales.  Therefore,  measurements  of  the  geometric  shape  of 
fracture  on  the  macroscopic  scale  can  be  used  to  infer  the  unit  process  on  the 
atomic  scale.  Likewise,  if  we  know  the  unit  process,  we  can  generate  a  shape. 
If  we  have  a  fractal  dimension  of  a  fracture  surface,  and  we  know  something 
about  the  structure  of  the  material,  then  we  can  deduce  possible  shapes  for  the 
unit  process.  We  could  then  build  a  fracture  surface  from  these  shapes. 


Our  second  assumption  involves  the  family  lines  in  Kj^  -  D*  space. 

These  lines  have  proposed  slopes  of  1/2.  The  slope  of  1/2  was  obtained 
from  observations  by  Mecholsky  et  al. 


This  assumes  that  all  of  the  materials  originally  tested  should  lie  along  the 
same  line  (i.e.,  are  members  of  the  same  family).  It  is  not  clear,  however,  that 
the  zinc  silicates  should  be  in  the  same  family  with  alumina. 

The  members  of  a  particular  family  should  exhibit  similar  fracture 
behavior.  Their  differences  in  toughness  should  reflect  differences  in  crack 
geometry,  not  fundamental  differences  in  fracture  mechanisms.  For  example, 
consider  a  hypothetical  material  which  is  toughened  by  the  addition  of  a 
second  phase.  This  material  would  have  three  basic  possibilities  for  crack  ' 
propagation,  only  one  bt  which  would  surely  place  it  in  the  same  family  as  the 
untoughened  material: 

(1)  Fracture  could  pass  through  the  second  phase  addition. 

(2)  Fracture  could  pass  through  the  interface  of  matrix  and 
second  phase., 

(3)  Fracture  could  deflect  around  the  second  phase  and  remain 
entirely  within  the  matrix. 

Only  the  third  possibility  would  represent  a  purely  geometric  toughening.  The 
crack  still  "sees"  the  same  material  yet,  because  of  deflection  around  the 
second  phase,  is  forced  to  follow  a  longer  path. 


A  series  of  heat  treated  Ocala  cherts  were  obtained  from  Ed 


Beauchamp  of  Sandia  National  Labs.  Purdy  and  Beauchamp^S  calculated  a 
decrease  in  toughness  as  a  result  of  heat  treatment.  Along  with  a  decrease  in 
toughness,  there  is  a  decrease  in  the  qualitative  fracture  surface  roughness. 

In  untreated  chert,  fracture  propagates  around  zones  of  high  density. 
Heat  treatment  improves  the  bonding  at  the  zone  boundariec  and  the  crack 
now  follows  a  less  tortuous  path  through  the  dense  zones.  It  may  seem 
paradoxical  that  improved  bonding  results  in  a  decrease  in  fracture  toughness; 
yet  the  improved  bonding  allows  the  crack  to  propagate  more  smoothly, 
thereby  generating  less  surface  area. 

The  dominance  of  geometrical  factors  in  affecting  the  toughness  of 
Ocata  chert  suggested  that  these  samples  would  constitute  a  family  in  the 

K(q-D*  plane.  Figure  27  shows  this  data  plotted  in  the  toughness-fractal 

dimension  plane.  A  least  squares  fit  gives  a  slope  that  is.  again,  0.5 .  Figure 
28  shows  coastline  sections  for  two  of  the  Ocala  chert  samples;  U5  and  #4. 
Note  the  dillerence  in  roughness  of  these  coastlines,  and  how  this  diflerence  is 
quanliiied  by  the  fractal  dimension. 


Another  assumption  that  went  into  the  determination  of  the  characteristic 
length  is  the  constant,  A,  of  equation  23.  The  components  of  A,  namely 
Young's  modulus  and  the  characteristic  length,  were  derived  from  a 
dimensional  analysis.  We  assured  the  existence  of  this  length  parameter  by 
assuming  that  it  was  reasonable  for  Young's  modulus  to  appear  in  equation 
(21).  This  is  hinged  upon  the  understanding  that,  ultimately,  for  a  material  to 
fracture,  bonds  must  break.  The  elastic  modulus  provides  an  average 
measure  of  the  strength  of  these  bonds  and.  so,  should  appear  in  a  model  of 


fracture. 


Pata._Summarv 


Two  techniques  were  used  to  determine  the  fractal  dimension: 

Slit-island  analysis  and  Richardson  plots.  Slit-island  analysis  uses 
area-perimeter  data  in  the  following  equation: 

A~p2yD 

Therefore,  logA  ~  2/D  logP. 

The  slope  of  a  graph  of  such  data  is  equal  to  2i'D,  where  A,  P.  and  D  are  a$ 
defined  earlier. 

Richardson  plots  employ  the  scale  dependence  of  measured  length  in: 

Therefoi^,  logL  =  (1  -D)  logE. 

The  slooe  of  such  a  relationship  is  equal  to  (t-D).  where  L  ,E.  and  D  are  as 
defined  earlier. 

Figures  29  through  41  contain  the  data  obtained  from  either  Richardson 
or  slit-island  analysis  for  the  materials  listed  in  Table  II.  ( This  excludes  the 
data  of  Wecholsky,  Feinberg  and  Passoja.) 
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Richardson  plot 
slope  =  1  -D 
D  =  1.15 


LogE 

Log  P 

2,72 

■n 

2.74 

3.06 

mmm 

3.16 

10.76 

3.53 

I  10.70 

Figure  31.  Data  for  F5  chert 


Richardson  plot 
slope  =  1-D 
D  =  1.26 
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Figure  32.  Data  for  chert  #3 


Figure  37.  Data  for  polycrystalline  spinel. 


Slit-island  plot 
slops  =  2/D 
D  =  1.16 
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Figure  38.  Data  lor  pyroceram 
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Richardson  plot 
slope  =  1-D 
D  =  1.20 
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Figure  39.  OalalorPZT. 


CHAPTER  4 
CONCLUSIONS 


The  irregular  structure  of  many  fracture  surfaces  presents  a  picture  that  is 
too  complicated  for  normal  geometric  description.  A  cross  section  through 
these  surfaces  will  provide  contour  lines  that  are  indicative  of  this  complicated 
geometry.  The  length  of  these  lines  is  related  to  the  roughness  and  area  of  the 
surface  from  which  it  was  derived.  Therefore,  classification  of  the  dimension  of 
such  a  line  is  useful  in  the  description  of  the  surface  from  which  it  was  obtained. 
The  basic  constructs  of  fractal  geometry  give  a  simple  methodology  for 
categorizing  the  tortuosity  of  a  line.  Thus,  fractals  are  used  to  classify  the 
geometry  of  a  line  passing  through  a  fracture  surface. 

A  relationship  is  shown  to  exist  between  a  material’s  fractal  dimension 
and  it«i  fracture  toughness; 


K|c  =  A(D-1)’'2. 


The  constant  A  is  as  a  family  parameter  that  groups  materials  in 

space  {where  0*  =  D*1 ).  The  value  of  this  constant  will  identify  a  particular 
family  line  that  may  represent  a  (racluro  mechanism. 
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Assuming  that  fracture  can  be  modeled  as  a  scaling  fractal,  we 
proposed  that  the  constant  A  is  a  product  of  Young's  modulus  and  a 

characteristic  length,  aQl 

A  =  E(ao)1'2 

Thus,  it  was  possible  to  compute  values  of  the  characteristic  length.  Modeling 
fracture  as  a  scaling  fractal  makes  it  possible  to  infer  the  atomic  and 
microscopic  geometric  structure  from  the  macroscopic  geometric  structure. 

Fractal  geometric  modeling  of  fracture  shows  a  number  of  encouraging 
results; 

(1)  The  dimension  ot  a  line  through  a  fracture  surface  gives  an 

indication  of  the  surface  toughness,  i.e.,  as  D  increases  the 
surface  roughness  increases. 

(2)  D  was  shown  to  be  related  to  fracture  toughness.  Within  a 
given  family:  as  D  increases,  toughness  increases. 

(3)  Modeling  the  fracture  surface  as  a  scaling  fractal  permits  the 
calculation  cl  a  characteristic  length  of  fracture. 


The  use  of  scaling  fractals  allows  for  an  orderly  construction  of  a 
complicated  geometric  shape.  Shapes  that  appear  too  chaotic  for  description 
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find  themselves  well  modeled  by  scaling  fractals.  Computer  simulated 
landscapes,  as  found  in  Mandelbrot’s  book,  are  strikingly  reminiscent  of  the 
natural  world.  It  seems  reasonable  that  these  fractal  models  will  be  useful  in 
describing  the  structure  of  the  objects  they  simulate. 

This  thesis  has  shown  that  a  fractal  dimension  can  be  assigned  to  brittle 
fracture  surfaces.  This  dimension,  as  yet,  is  restricted  to  the  description  of  ‘ 
contour  lines  on  the  fracture  surface,  as  well  as  fracture  sudace  profiles.  It  must 
be  extended  to  the  actual  two-dimensional  nature  of  the  surface.  The  surface 
dimension  will  lie  between  two  and  three-two  being  a  penectly  plane  fracture. 

As  the  dimension  increases  from  two  to  three,  vhe  surface  will  be  increasing  in 
roughness  and  geometiic  complexity.  The  relation  between  fractal  dimension, 
fracture  toughness,  fracture  surface  roughness  and  suriace  area  cannot  be 
properly  explored  without  the  two-dimensional  couriterpart  of  the  dimensions 
determined  in  this  thesis. 

It  is  clear  that  fracture  models  have  taken  no  account  of  the  tremendous 
complexity  ol  fracture.  In  fac*,  the  moaels  assume  that  a  crack  simply  separates 
planes  of  atoms  into  two  perfectSy  plane  fracture  surfaces.  The  fractal  approach 
will  enable  the  accommodation  of  surface  roughness  into  the  basic  models. 
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Instead  of  separating  planes  of  atoms,  the  crack  is  aliowed  to  wander. 
The  new  profiles  generated  by  such  a  crack  will  be  composed  of  shapes  that 
resemble  fractal  generators.  Thus,  fracture  surfaces  can  be  constructed  from 
atomic  scale  shapes  that  cascade,  via  the  scale  invariance  of  fractal  geometry, 
'■j  macroscopic  fracture  surface  features.  Computer  simulations  of  this 
hypothesis  are  the  next  step  in  understanding  the  applicability  of  fractal 
concepts  to  brittle  fracture. 
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